Parity violating superfluidity in ultra-cold fermions 
under the influence of artificial non-Abelian gauge fields 
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We discuss the creation of parity violating Fermi superfluids in the presence of non-Abelian gauge 
fields involving spin-orbit coupling and crossed Zeeman fields. We focus on spin-orbit coupling 
with equal Rashba and Dresselhaus (ERD) strengths which has been realized experimentally in 
ultra-cold atoms, but we also discuss the case of arbitrary mixing of Rashba and Dresselhaus (RD) 
and of Rashba-only (RO) spin-orbit coupling. To illustrate the emergence of parity violation in the 
superfluid, we analyze first the excitation spectrum in the normal state and show that the generalized 
helicity bands do not have inversion symmetry in momentum space when crossed Zeeman fields are 
present. This is also reflected in the superfluid phase, where the order parameter tensor in the 
generalized helicity basis violates parity. However, the pairing fields in singlet and triplet channels 
of the generalized helicity basis are still parity even and odd, respectively. Parity violation is 
further reflected on ground state properties such as the spin-resolved momentum distribution, and 
in excitation properties such as the spin-dependent spectral function and density of states. 

PACS numbers: 03.75.Ss, 67.85.Lm, 67.85.-d 



Parity violating phenomena are very rare in physics, 
but a classical example is known from particle physics, 
where parity violating processes of the weak interaction 
were proposed pj and observed in the decay of 60 Co sev- 
eral decades ago 0. In this case, the weak interactions 
allow for parity violation, but the particle kinetic energies 
are parity even, reflecting the inversion symmetry of their 
space. The Standard Model of particle physics, which is 
a non-Abelian gauge theory, incorporates parity viola- 
tions and postulates that for nuclear beta decay parity is 
maximally violated. Other examples of parity violation 
exist for instance in condensed matter physics, where par- 
ity breaking is associated with crystals without inversion 
symmetry [3[ or with crystals which have inversion sym- 
metry initially, but can develop spontaneously permanent 
electric polarization through lattice distortions leading 
to ferroelectric materials [^|. However, examples of par- 
ity breaking in superfluids, such as those encountered in 
nuclear, atomic, condensed matter and astrophysics are 
hard to find, and to our knowledge there seems to be no 
confirmed example in nature. 

Recently, it has been possible to create non-Abelian 
gauge fields in ultra-cold atoms via artificial spin-orbit 
(SO) coupling of equal superposition of Rashba [j| 
hfl(k) = Vr(— k y x + k x y) and Dresselhaus @ hrj(k) = 
voikyX. + k x y) terms, leading to the equal- Rashba- 
Dresselhaus (ERD) form 0, [|| 1o.erd0<-) = vk x y, where 
vr = vd = for which parity preserving superfluid- 
ity is possible |9HL1|. Other forms of SO fields, such as 
the Rashba-only or Dresselhaus-only cases, require ad- 
ditional lasers and create further experimental difficul- 
ties [l2|, while several theory groups have investigated 
the Rashba-only case Hl-liJ due to the connection to 
earlier condensed matter literature 17-19^ 

The current Zeeman-SO Hamiltonian created in the 



laboratory is 
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zso 



(k) 



[hy + hERDQt)]cr y (1) 



for an atom with center-of-mass momentum k and spin 
basis | f), | I). The fields h z = -Qr/2, h y = -8/2, and 
^■Bijn(k) = vk x can be controlled independently. Here, 
Q r is the Raman coupling and S is the detuning, which 
can be adjusted to explore phase diagrams as achieved in 
87 Rb experiments Q, or to study the high-temperature 
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normal phases of Fermi atoms [2C 

In this letter, we show that ultra-cold Fermi superflu- 
ids in the presence of non-Abelian gauge fields consist- 
ing of artificially created spin-orbit and crossed Zeeman 
fields described in Eq. ([T]) can produce a parity violating 
superfluid state when interactions are included. How- 
ever, unlike the case of the Standard Model where parity 
breaking is driven by the weak force, in our case, par- 
ity breaking is driven by the effects of the non-Abelian 
gauge field on the kinetic energy. To illustrate the lack of 
parity in physical observables, we analyze spectroscopic 
quantities such as the elementary excitation spectrum, 
momentum distribution, spectral function and density of 
states in the superfluid state. 

Hamiltonian: To analyze parity violation in ultra-cold 
Fermi superfluids, we start from the Hamiltonian in mo- 
mentum space as 



H n 



5>J(k)Ho(kty.(k), 



(2) 



ks 



where H (k) = [K(k)l - h e s(k) • a] with K(k) = 
k 2 /2m — /i being the single particle kinetic energy rela- 
tive to the chemical potential fi; the vector-matrix a de- 
scribes the Pauli matrices (a x , cr y , a z ); h e g (k) is the effec- 
tive magnetic field with components [h x (k), h y (k), h z (k)} 
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and V4(k) is the creation operator for fermions with 
spin s and momentum k. In the ERD case, which is 
readily available in ultra-cold atoms, the effective mag- 
netic field is simply h cff (k) = [0 7 h y + h ERD (k),h z ] 7 
where h y and h z are Zeeman components correspond- 
ing to the detuning 5 and the Raman coupling fi/j, while 
^EiiD(k) = vk x is the spin-orbit field. We define the 
total number of fermions as N = iV-j- + iVj., and the in- 
duced population imbalance as P; n d = (JVf — Ni)/N. Wc 
choose our scales through the Fermi momentum fc^ de- 
fined from N/V = k F /(3ir 2 ), leading to the Fermi energy 
(F = k E /2m and the Fermi velocity Uj? = kp/m. 

Generalized Helicity Basis: The matrix Ho(k) can 
be diagonalized in the generalized helicity (GH) basis 
k, a) = $t (k)|0) via a momentum-dependent SU(2) ro- 
tation generated by the unitary matrix 



U k = 



-«k w k 



(3) 



where the normalization condition |uk| 2 + |^k| 2 = 1 is 
imposed to satisfy the unitarity condition U^Uk = 1. 
The corresponding eigenvectors are the spinors $(k) = 
Uj ( .*(k), where $(k) = [^(k), "^(k)] is expressed in 
terms of ip(k) — [V't(k), ^(k)] by the relations $f|-(k) = 
^kCkt — ^kCkj. and $j).(k) = v£c k t + itkCk|- The coherence 



factor ttk 
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|h off (k) 



is chosen to be real with- 



out loss of generality, and «k = —e l 



1 - 



h of f(k) 



is a complex function with phase defined by </?k = 
Arg[/ij_(k)] . The complex field h±(k) = h x (k) — ih y (k) 
has components h x (k) and h y (k) along the .t and y direc- 
tions, respectively. The magnitude of the effective field is 
|h eff (k)| = yjh 2 z + |Mk)| 2 . In the ERD case h x (k) = 0, 
and the ratio hx(k)/\h x (k)\ = e lVk = -isgn \h y (k)} , 
where h y (k) = h y + vk x . 

The generalized helicity spins a = JJ-) are aligned 
or antialigned with respect to the effective magnetic field 
h e ff(k), and the corresponding eigenvalues of H (k) are 
£fr(k) = e^(k) - fx and £^(k) = ejj.(k) - \i. Here, the he- 
licity energies are simply e^-(k) = K(k) — \h c g(k)\ and 
ej).(k) = K(k) + |/i c ff(k)|. In the specific case of ERD 
coupling with non-zero detuning (h y ^ 0) the effective 
field is h c ff (k) = h z i + [h y + IierdO<)} y, with magni- 
tude |/i c g (k)| = \J h? z + (h y + vk x ) 2 and parity violation 
occurring along the x axis. This is illustrated in Fig. [TJ 
where for finite h y (non-zero detuning 5) the generalized 
helicity bands e^(k) and e^(k) do not have well defined 
parity in momentum space. As seen in Fig.[lja)-(b), par- 
ity is preserved for v ^ if h y = (zero detuning). While 
as noted in Fig. HKc)-(d), parity is violated for v ^ 0, if 
h y (finite detuning). Similar parity violation along 
the x axis occurs for other mixtures of Rashba and Drcs- 
sclhaus terms as long as h y ^ 0. 
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FIG. 1. (color online) Generalized helicity bands e^(k)/eF 
(blue line) and e^(k)/eF (red line) versus momentum k x /kp 
with k y — k z — and for ERD spin-orbit coupling v/vf = 
0.4. The black dashed lines show the helicity bands for 
v/vf = 0.4 with h z /£F = h y /eF = 0. The Zeeman fields 
are (a) h y /eF — and h z /eF = 0.1, (b) h y /eF = and 
h z /eF = 0.7, (c) h y /eF = 0.1 and h z /eF = 0.1, (d) h y /eF = 
0.2 and h z /eF = 0.7. Notice that e Q (k) ^ e Q (— k) in (c) and 
(d), indicating the absence of parity. 



Interactions and Order Parameter: In order to under- 
stand the underlying physics of this system, it is im- 
portant to rewrite the interaction Hamiltonian in the 
generalized helicity basis. The starting interaction is 
rii = -sEq^tl)^)! where the pair creation op- 
erator with center of mass momentum q is 6^(q) = 
X^k ■^l(k+q/2)'0| (— k+q/2), can be written in the helic- 
ity basis as Hi = -g^qap-yS -^(qJ-frysfa)' where the 
indices a,/3,7,<5 cover fr and JJ. states. Pairing is now 
described by the operator 



-B Q/ 9(q) 



k 



A^(k + ,k_)$ Q (k+)^(k_) (4) 



and its Hermitian conjugate, with momentum indices 
k± = ±k + q/2. The matrix A Q( g(k + , k_) is directly re- 
lated to products of coherence factors u(k±), v(k±) (and 
their complex conjugates) of the momentum dependent 
SU(2) rotation matrix U(k±). Seen in the GH basis, the 
interactions reveal that the center of mass momentum 
k + + k_ = q and the relative momentum k + k = 2k 
are coupled and no longer independent, and thus do 
not obey Galilean invariance. The interaction constant 
g is related to the scattering length via the Lippman- 
Schwinger relation V/g = —Vm/(4na s ) + J2k V(2 £ k)- 

From Eq. (j4]) it is clear that pairing between fermions 
of momenta k+ and k_ can occur within the same helic- 
ity band (intra-helicity pairing) or between two different 
helicity bands (inter- helicity pairing) . For pairing at zero 
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center-of-mass momentum q = 0, the order parameter 
for superfluidity is the tensor A a p(\t) = A n A a p(k, — k), 
where A = —g J ^2 jS (B 1 s(0)), leading to components: 
A^(k) = A (wk^-k — «k u -k) f° r total helicity pro- 
jection A = +1; A^(k) = -A (u k u-k + Wk"l k ) and 
Ajj,^(k) = Ao (ukW-k + ^k^-k) f° r total helicity projec- 
tion A = 0; and A,y.|j.(k) = Ao (u^v^ — v^u-k) for total 
helicity projection A = — 1. Parity is violated in A Q , ( g(k) 
since they do not have well defined parity for non-zero 
spin-orbit coupling and crossed Zeeman fields h y and h z . 

However, we may still define singlet and triplet sec- 
tors in the generalized helicity basis, which are even and 
odd in momentum space respectively for any value of h y . 
The singlet sector is defined by the scalar order param- 
eter A5, (k) = [A w (k)-A w (k)] /2 corresponding to 
A = 0. While the triplet sector is defined by the vector 
order parameter A^^k), by its generalized helicity com- 
ponents A<r,+i(k) = A^(k) corresponding to A = +f; 
A T ,o(k) = [A^(k) + A^(k)] /2 corresponding to A = 0; 
AT,-i(k) = A^(k) corresponding to A = — f. 

Superfluid Ground State and Elementary Excita- 
tions: The ground state for uniform superfluidity 
can be expressed in terms of fcrmion pairs in the 
GH basis as the many-body wavefunction \G) = 

U a p(k) + ^(k)$t(k)*J(-k)] } |0), where 

0) is the vacuum state with no particles. 
The Hamiltonian matrix in the GH basis is 



n k {E c 



H cx (k) 



/ ^ A n (k) A u (k)\ 

^ A M (k) %(k) 

A* u (k) A* u (k) -£_ kt 

\A* n (k) A*^(k) -e-^/ 



which is traceless, showing that the sum of its eigenval- 
ues is zero. We have obtained analytical solutions for 
the eigenvalues of H cx (k) for arbitrary RD spin-orbit or- 
bit and arbitrary Zeeman fields h y and h z , but we do 
not list them here, because their expressions are cumber- 
some. However, for each momentum k, the determinant 

Det tul — H cx (k) , leads to the quartic equation 

w 4 + a 3 (k)w 3 + a 2 (k)w 2 + ai(k)w + a (k) = 0. (6) 

In the particular case of ERD spin-orbit coupling with 
crossed Zeeman fields, the coefficients become 03 (k) = 0, 
the coefficient of the quadratic term takes the form 



o 2 (k) = -2 (K 2 (k) + |A | 



\vk r , 



\hyV + \K 



while the coefficient of the linear term is ai(k) = 
— 8K(k)(vk x )h y , and lastly the coefficient of the zero-th 
order term is 

oo(k) = ^(k)^(k)^(-k)^(-k) + |A | 2 a 2 (k), 

where a 2 (k) = (2K 2 (k) + |A | 2 + h 2 (k)) with h 2 Q (k) = 
2\vk x \ 2 - 2\h y \ 2 - 2\h z \ 2 . Notice that a 2 (k) and a (k) 



have even parity, while a\(k) has odd parity and is thus 
responsible for the parity violation that occurs in the ele- 
mentary excitation spectrum. Furthermore, parity viola- 
tion occurs only when both v and h y are non-zero, since 
when either h y = or v = the coefficient a\ (k) vanishes 
and parity in the elementary excitation spectrum is fully 
restored. From the secular equation, it follows that when 
k x = 0, the coefficient ai(k) also vanishes and the excita- 
tion energies Ei(0, k y ,k z ) have the same analytical form 
as in the case for h y = 0, with the simple replacement 



of hi — > h' 2 + h 2 . This property is just a consequence 
of the reflection symmetry of the Hamiltonian through 
the k x = plane. However, parity is violated, be- 
cause inversion symmetry through the origin of momenta 
does not exist, that is, Ei{— k) ^ Ei(k). In contrast, 
quasiparticle-quasihole symmetry is preserved since the 
corresponding quasiparticle-quasihole energies obey the 
relations E 2 {k) = -E 3 (-k) and £d(k) = -E 4 (-k). 

A simple inspection shows that gapless and fully 
gapped phases emerge. A gapless phase with two rings 
of nodes (US-2) appears when h 2 + h 2 — |A | 2 > 

and /i > 4 h 2 + h% — |Aq| 2 . A gapless phase with one 



ring of nodes (US- 1) occurs 
and \n 



for hy 



hi - I An 



> 



< 



h 2 — |A | 2 . A directly gapped phase 



(d-US-0) arises for hi 



K - I Ar 



> and 



h? y 



h 2 — I Ao 1 2 , while an indirectly gapped phase (i- 
US-0) emerges for hi + h 2 z - | A | 2 < and [i > 0. Lastly, 



(5) the quasiparticle excitation energy -E 2 (k) becomes neg- 



ative in certain momentum regions when h y > |Ar 
indicating that the uniform ground state becomes less 
energetically favorable against the normal state 22 1 . 

Phase Diagram and Thermodynamic Poten- 
tial: From the thermodynamic potential fl\js = 
-(7 1 /2)Ek,,fo[l + cxp(-^(k)/T)] + Ek^(k) + 
Ao| 2 /g we obtain self-consistently the zero tempera- 
ture (T = 0) phase diagram as a function of crossed 
Zeeman fields h y and h z for vjvp = 0.4 at unitarity 
l/(kFd s ) = in Fig. [U[a), and at the BEC regime 
f/(fci?a s ) = 2.0 in Fig. [2jb), but a stability analysis 
against non-uniform phases is necessary as in the 
parity- preserving case [ToL \vj\ . At unitarity the uniform 
superfluid phases i-US-0, US-1, US-2 and the normal 
(N) phase are present in the range shown, while in the 
BEC regime only the d-US-0 occurs in the same range 
of fields. The transitions between different US phases 
is topological with no change in symmetry as in the 
parity- preserving case 11|- While the transitions 
from US phases to the N phase involve a change in 
symmetry, from broken to non-broken U(I), and are 
discontinuous, as seen in the insets of Fig. [2j 

Detecting parity violation: A direct measurement 
of parity violation in the superfluid state can be 
made through the momentum distributions n s (k) = 
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FIG. 2. The T = phase diagram in the h y -h z parameter 
space showing various uniform superfluid phases US-2, US-1, 
d-US-0 and i-US-0, and the normal phase for ERD spin-orbit 
coupling v/vf — 0.4 and at (a) unitarity l/(fcj^a s ) = 0.0 
and in (b) the BEC regime l/(fcj^a s ) = 2.0. The insets show 
|Ao| as a function of h y for h z /eF = 0.2 (dotted line); for 
hz/tF = 0.4 (dot-dashed line); h z /eF = 0.6 (dashed line); 
and hz/eF = 0.8 (solid line). In the range shown, |Ao| is 
essentially independent of h y and h z in the BEC regime. 



(V4(k)iMk)). They arc illustrated in Fig. [3] for US- 
1 superfluid with spin-orbit v/vp — 0.4 and interac- 
tion l/(kpa s ) = 0, in the parity-preserving case with 
hy/eF — and h z /ep = 0.7 in (a)-(d) and in the parity- 
violating case with h y /ep = 0.2 and h z /ep = 0.7 in (e)- 
(h). At finite temperatures, the momentum distributions 
broaden, but parity violation is still self-evident. 
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FIG. 3. (color online) Momentum distributions (T = 0) 
w-f(k) (two left-most columns) and nj.(k) (two right-most 
columns) for l/(kFCis) = 0.0 and v/vf = 0.40 at the US-1 
phase. In (a)-(d) the field values are h y /eF = 0, h z /eF = 0.7, 
with (i/e F = 0.5803, |A |/e F = 0.3592, and P ind = 0.6592. In 
(e)-(h) the field values are h y /eF = 0.2 and h z /eF = 0.7, with 
p/e F = 0.5871, |A |/e F = 0.3157, and P ind = 0.6958. The 
blue-dashed and red-solid lines represent cuts of n 3 (k) along 
the directions (0,k y ,0) and (k x ,0,0), respectively. 

Parity violation is also manifested in other momen- 
tum resolved properties such as the spectral func- 
tion A s (u>,)s.) = — (l/Tr)lmG ss (iuj = u> + iS, k), where 



G ss (iuj, k) = iuil — -ffex(k) 



written in the 



--t,i 



festation of parity violation in the elementary excitation 
spectrum for the US-1 superfluid phase, and the corre- 
sponding implications for momentum integrated quanti- 
ties such as the spin- resolved density of states p s {uS) = 
A s (io, k). The most important point is that for finite 
spin-orbit coupling v and when h y ^ 0, the excitation 
energies Sj(k) ^ Ei(—\c). This implies that degenerate 
peaks at h y = (corresponding to minima or maxima of 
the excitation spectrum) are increasingly split with grow- 
ing h y . This effect is illustrated in Fig. @] at the locations 
indicated by the small black arrows. 
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basis. Instead, in Fig. [H we choose to illustrate a mani- 



FIG. 4. (color online) Eigenvalues Ei(k) and density of states 
ps{u) (in units of tF and e F , respectively) for l/(fci?ffls) = 
and v/vf = 0.4 in the US-1 phase, but close to the US-l/US-2 
boundary, with parameters h y /eF = 0, h z /eF = 0.7, p/eF = 
0.5803, |A |/e F = 0.3592, and P ind = 0.6592 in (a)-(d); and 
hy/ep = 0.2, hz/ep = 0.7, p/e F = 0.5871, |A |/eF = 0.3157, 
and P illd = 0.6958 in (e)-(h). Cuts of E % (k) along (fee, 0,0) 
are shown in (a) and (e), and along (0, k y , 0) are shown in (d) 
and (h). In panels for p s (oj) a small broadening 8/ef = 0.01 
is used. The black arrows indicate examples of peaks that 
split when parity breaking occurs for finite h y . 



Conclusions: We showed that non-Abelian gauge fields 
consisting of spin-orbit and crossed Zeeman fields lead 
to parity violating superfluidity in ultra-cold Fermi sys- 
tems. We derived general relations that can be applied 
to spin-orbit couplings involving any linear combination 
of Rashba and Dresselhaus terms. We focused mostly on 
the case of equal Rashba-Dresselhaus (ERD) spin-orbit 
coupling. The presence of such fields produce a super- 
fluid order parameter tensor whose components in the 
generalized helicity basis are neither even nor odd un- 
der spatial inversion. Even though the elements of this 
tensor written in generalized singlet or triplet helicity 
channels have even or odd parity, respectively, the exci- 
tation spectrum does not have well defined parity, but 
preserves quasiparticle-quasihole symmetry. This parity 
violation has important experimental signatures leading 
to momentum distributions without inversion symmetry 
and to spin-resolved density of states that possess split 
peaks in frequency. 
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